% function [baseMVA, bus, gen, branch, genpara, NTswitch, sbspara,Hess_coef,u0,mu_coef] = case9
baseMVA = 100.00;

% bus  type1 type2  Pdemand  Qdemand   Gshunt   Bshunt   Vmodule    Vangle  baseKV  Vmax  Vmin
% type1: PQ bus = 1; PVbus = 2; reference bus = 3; isolated bus >= 4
% type2: non-load = -1; load -- no shedding = 0; load -- shedding =1
bus = [
     1     3    -1      0.0000      0.0000      0.0000      0.0000      1.0400      0.0000    345.0000      1.1000      0.9000
     2     2    -1      0.0000      0.0000      0.0000      0.0000      1.0250      9.2800    345.0000      1.1000      0.9000
     3     2    -1      0.0000      0.0000      0.0000      0.0000      1.0250      4.6650    345.0000      1.1000      0.9000
     4     1    -1      0.0000      0.0000      0.0000      0.0000      1.0258     -2.2170    345.0000      1.1000      0.9000
     5     1     1    125.0000     50.0000      0.0000      0.0000      0.9956     -3.9890    345.0000      1.1000      0.9000
     6     1     1     90.0000     30.0000      0.0000      0.0000      1.0127     -3.6870    345.0000      1.1000      0.9000
     7     1    -1      0.0000      0.0000      0.0000      0.0000      1.0258      3.7200    345.0000      1.1000      0.9000
     8     1     1    100.0000     35.0000      0.0000      0.0000      1.0159      0.7280    345.0000      1.1000      0.9000
     9     1    -1      0.0000      0.0000      0.0000      0.0000      1.0324      1.9670    345.0000      1.1000      0.9000
];

% bus       Pg         Qg     Vsetpoint  base   status1 status2[cuttable?] Pmax  Pmin  Qmax  Qmin
gen = [
     1     71.6410     27.0500      1.0400    100.0000     1     1    250.0000     10.0000    300.0000   -300.0000
     2    163.0000      6.6500      1.0250    100.0000     1     1    300.0000     10.0000    300.0000   -300.0000
     3     85.0000    -10.8600      1.0250    100.0000     1     1    270.0000     10.0000    300.0000   -300.0000
];

% fbus  tbus  r   x    b      ratea    rateb     ratec Vf/Vt-ratio  angle status[offline?]
% rateA: long term rating; rateB: short term rating; rateC: emergency rating
branch = [
     1     4      0.0000      0.0576      0.0000    250.0000    250.0000    250.0000      1.0000      0.0000     1
     4     6      0.0170      0.0920      0.1580    250.0000    250.0000    250.0000      0.0000      0.0000     1
     6     9      0.0390      0.1700      0.3580    150.0000    150.0000    150.0000      0.0000      0.0000     1
     3     9      0.0000      0.0586      0.0000    300.0000    300.0000    300.0000      1.0000      0.0000     1
     8     9      0.0119      0.1008      0.2090    150.0000    150.0000    150.0000      0.0000      0.0000     1
     7     8      0.0085      0.0720      0.1490    250.0000    250.0000    250.0000      0.0000      0.0000     1
     7     2      0.0000      0.0625      0.0000    250.0000    250.0000    250.0000      1.0000      0.0000     1
     5     7      0.0320      0.1610      0.3060    250.0000    250.0000    250.0000      0.0000      0.0000     1
     5     4      0.0100      0.0850      0.1760    250.0000    250.0000    250.0000      0.0000      0.0000     1
];

% bus   model  startup  shutdown n_coef   Xd'       H=1/2*Tj  D
% model: 1 -- piecewise linear   2 -- polynomial
genpara = [
     1     2   1500.0000      0.0000     3      0.0610     23.6400      0.0000
     2     2   2000.0000      0.0000     3      0.1200      6.4000      0.0000
     3     2   3000.0000      0.0000     3      0.1810      3.0100      0.0000
];

% contingency set
% Type    Loctn1   Loctn2    Occurr      Clear_fault    Clear_load     R           X
NTswitch = [
    1110     7     5      0.0000      0.1700      0.3000  1.0000e-010  1.0000e-010
  %1110     7     5      0.0000      0.1700      0.3000  1.0000e-010  1.0000e-010
  %1110     7     8      0.0000      0.2000      0.3000  1.0000e-010  1.0000e-010
];

% parameters for transient stability
sbspara = [
        0.01           % step size
        100 *pi/180    % radian of max angle offset
       60.00           % frequency
];
